Signatures of valence fluctuations in CeCu 2 Si2 under high pressure 
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Simultaneous resistivity and a.c.-specific heat measurements have been performed under pressure 
on single crystalline CeCu2Si2 to over 6 GPa in a hydrostatic helium pressure medium. A series of 
anomalies were observed around the pressure coinciding with a maximum in the superconducting 
critical temperature, T™ ax . These anomalies can be linked with an abrupt change of the Ce valence, 
and suggest a second quantum critical point at a pressure P v ~ 4.5 GPa, where critical valence 
fluctuations provide the superconducting pairing mechanism, as opposed to spin fluctuations at 
ambient pressure. Such a valence instability, and associated superconductivity, is predicted by an 
extended Anderson lattice model with Coulomb repulsion between the conduction and f-electrons. 
We explain the T-linear resistivity found at P v in this picture, while other anomalies found around 
P v can be qualitatively understood using the same model. 



I. INTRODUCTION 

Boosted by the discovery of superconductivity in 
CeCu2Si2 over twenty years agoji the relationship be- 
tween superconductivity and magnetism has been ex- 
tensively investigated in various d and f heavy fermion 
(HF) compounds. A consensus has developed that HF 
superconductivity is mediated by spin fluctuations; 2 ' 3 ' 4 ' 5 
mainly because superconductivity was found close to a 
magnetic instability at T = 0, sometimes described as a 
quantum-critical point (QCP), often attained by apply- 
ing pressure. A recent development is that an essentially 
gapless superconducting (SC) state has been identified 
by NMR/NQR measurements in the region where the 
SC state coexists with antiferromagnetismj&l consistent 
with a theoretical prediction. 8 

In this paper we further explore the possibility that 
at a pressure P v ~ 4.5 GPa, a second QCP, associated 
with f-electron occupation number, has a major role to 
play in the superconductivity of CeCu2Si2 and related 
compounds Ai£ 

CeCu2Si2 has a superconducting ground state at ambi- 
ent pressure with a critical temperature T c , around 0.7 K. 
It is firmly believed that the compound is close to an an- 
tiferromagnetic QCP at slight negative pressure, accessi- 
ble for example by partial substitution of Si with Gc. 11 
When pressure is applied, T c initially remains close to its 
ambient pressure value, followed by a sudden increase to 
around 2 K at about 3 GPa. Further increase in pressure 
results in a slower suppression of T c to zero. This non- 
monotonous behavior of T C {P) was first explored by re- 
sistivity in the quasi-hydrostatic conditions of Bridgman 
anvil cell. 12 Subsequent investigations by susceptibility 13 
and resistivityi 4 * 1 ^ were carried out in various pressure 
media, and showed considerable variation in T c between 
samples, especially at high pressure. 

With increasing pressure CeCu2Si2 passes from a 



nearly trivalent 4f 1 behavior, with Kondo coupling be- 
tween conduction and f-electrons, to behavior at very 
high pressure characteristic of intermediate valence (IV) 
systems, whose valence fluctuates between the 4f™ and 
4f n ~ 1 + [5d6s] electronic configurations. As a result, deep 
in this IV regime, the resistivity, for instance, resembles 
that of LaCu2Si2, which lacks 4f electrons. Similar T C (P) 
dependence to that found in CeCu2Si2 is seen in the iso- 
electronic sister compound, CeCu2Ge2, offset by about 
10 GPa due to the larger atomic volume of Ge>i& Apart 
from this shift of the pressure scale, no qualitative dif- 
ference occurs between the two compounds beyond un- 
avoidable sample dependence. 

From a more theoretical point of view, there exist at 
least three reasons to believe that critical valence fluctu- 
ations are at the origin of the pressure-induced peak of 
the SC transition temperature T c . 

First, the A coefficient of the T 2 resistivity law de- 
creases drastically by about two orders of magnitude 
around the pressure corresponding to the T c peak. 10 
Since A scales as (m* /m) 2 in the so-called Kondo regime, 
this implies that the effective mass m* of the quasiparti- 
cles also decreases sharply there. This fall of m* is possi- 
ble only if there is a sharp change of Ce valence, deviating 
from Ce 3+ , since the following approximate formula for 
the renormalization factor q holds in the strongly corre- 
lated limit iiii 8 . 



m* 1 1 — rtf/2 

m 1 — nf 



(1) 



where ri{ is the f-electron number per Ce ion. 

Second, the so-called Kadowaki- Woods (KW) ratio^ 
A/7 2 , where 7 is the Sommerfeld coefficient of the elec- 
tronic specific heat, crosses over quickly from that of a 
strongly correlated class to a weakly correlated oneiSi 
The inverse of the Sommerfeld coefficient, 7 _1 , can be 
identified with the Kondo temperature Tk, which is 
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experimentally accessible by resistivity measurements. 
This indicates that the mass enhancement due to the 
dynamical electron correlation is quickly lost at around 
P ~ P v , in agreement with the previous point. The phe- 
nomenon can be understood if we note the fact that 7 
consists essentially of two terms: 

7 = 7band I 1 7^— I ) 

= 7band+7cor, (2) 

where 7band is due to the so-called band effect and 7 cor = 
— 7band9E(e)/i9e is due to the many-body correlation ef- 
fect, with E(e) being the self-energy of the correlated 
electrons. j cot and A are related to each other through 
the Kramers-Kronig relation, leading to the large value 
of the KW rations and when 7 cor ^> 7band, this is in- 
deed seen. On the other hand, if 7 CO r ~ 7band, the ratio 
A/7 2 should be reduced from the KW value consider- 
ably because the effect of 7band cannot be neglected in 
its denominator. 

Third, there is a sharp peak in the residual resistiv- 
ity po at around P ~ P v r^ which can be understood 
as a many-body effect enhancing the impurity potential 
(in fact we define the pressure P v experimentally by the 
maximum of po). In the forward scattering limit, this 
enhancement is proportional to the valence susceptibility 
— (drif/def)^, where C{ is the atomic f- level of the Ce ion, 
and p is the chemical potential, 21 Physically speaking, 
local valence change coupled to the impurity or disorder 
gives rise to a change of valence in a wide region around 
the impurity which then scatters the quasiparticles quite 
strongly, leading to the increase of pq. The enhancement 
of po can be thus directly related to the degree of sharp- 
ness of the valence change, because the variation of the 
atomic level ef is considered to be a smooth function of 
the pressure. 

These circumstantial clues to the importance of crit- 
ical valence fluctuations have been backed up by a mi- 
croscopic calculation of T c for d- wave pairing as a func- 
tion of ef iSk This showed that sudden valence change oc- 
curs if a moderately sized Coulomb repulsion U c f is taken 
into account between the conduction c-, and localized f- 
electrons, with the peak structure of T c being qualita- 
tively reproduced. 

Table [I] summarizes the current experimental evidence 
of anomalies seen in CeCu2(Ge/Si)2 around P v . 

Part (i) of table HI refers to direct evidence for a va- 
lence transition of the Ce ion: Cell volume^ 3 , and Ljjj 
X-ray absorption 24 measurements show discontinuities as 
a function of pressure. The drastic decrease of the A co- 
efficient of the T 2 resistivity law, along with the A vs 
j>max sca ii n g relation, indicate that the system is leav- 
ing the strongly correlated regime characterized by a f- 
occupation number close to unity. (T x max is defined in 
Fig. Eland assumed to be proportional to Tr-.) 

Part (ii) refers to anomalies observed close to the 
maximum of T c predicted by critical valence fluctuation 



theory. 21,2 — These are the maximum of T c itself and the 
enhanced residual resistivity, p J^i^i^i^i^i^ 

Part (iii) refers to properties following from the ex- 
tended treatment of the critical valence fluctuations 
found in section lTVI of this paper. This includes the linear 
resistivity^ and the maximum in 7, both found around 
P.- 

In part (iv) are listed the remaining features that are 
observed in CeCu2Si2 and CeCu2Ge2 around the maxi- 
mum in T c but which are so far not fully explained. For 
example the merging of T 1 nax and T 2 max , where the latter 
(also defined in Fig. is believed to reflect the effect of 
the excited crystalline electric field (CEF) split f-levels. 
Many of the anomalies noted in table^have also been ob- 
served to coincide with the maximum of T c in other HF 
superconductors, from CePd 2 Si2 29 to CeCusAu^, the 
latter showing traces of superconductivity under pres- 
sure. 

Previous work on CeCu2Si2 has shown a lot of varia- 
tion in low-temperature behavior between different sam- 
ples. The fact that not all reports have shown ev- 
ery anomaly is not entirely surprising, since large vari- 
ations in the electronic properties of CeCu2Si2 are well 
known to result from extremely small differences in 
compositioni^ 1 *^ The extension of these variations with 
pressure has not been systematically explored, but al- 
most all samples so far studied have shown an enhance- 
ment of T c , along with effects such as the enhancement 
of the residual resistivity, to be discussed below. 

This variability under pressure may be due to the sam- 
ples themselves, or to pressure inhomogeneities caused by 
non-hydrostatic pressure media. We were therefore mo- 
tivated to use solid helium as a pressure medium, due to 
its near-ideal hydrostaticity at low temperature. By si- 
multaneously probing resistivity and specific heat in the 
same sample, we were able to explore both percolative 
transport and bulk evidence for superconductivity. 



II. EXPERIMENTAL METHODS 

High pressure was induced using diamond anvils with 
a 1.5 mm culeti^ A stainless steel gasket was specially 
prepared to absorb the large volume decrease of the he- 
lium pressure medium from ambient pressure, and to 
avoid severing the measurement wires. These were insu- 
lated from the gasket using a mixture of AI2O3 powder 
and epoxy resin. The pressure was measured to within 
0.02 GPa at various temperatures down to 4.2 K using 
the ruby fluorescence scale. 

The CeCu2Si2 sample was prepared by reaction of its 
constituent elements with a slight excess of Cu, with a 
nominal initial composition CeCu2.iSi2. The product 
was then melted in an induction furnace and slowly al- 
lowed to crystallize under 50 bars Ar in a BaZrOs crucible 
(see Ref. [lj for more details). 

The small monocrystal used in this work was cut and 
polished to 230 x 80 x 20 /im 3 , and six 5 pm <f> wires (four 
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TABLE I: Anomalies in CeCu 2 Si 2 and CeCu 2 Ge 2 associated with valence transition, with references. Symbols explained in 
the text. 

Part (i): Direct evidence for sudden valence change. 

Part (ii): Anomalies explained by published valence fluctuation theoryi— 

Part (iii): Anomalies explained by extended treatment of the critical valence fluctuations (section llVI . 
Part (iv): Other anomalies observed around crossover to intermediate valence with pressure. 



gold and two Au + 0.07% at. Fe) were spot welded to 
the sample. The c-axis of the tetragonal structure was 
parallel to its smallest dimension. The magnetic field, 
when applied, was parallel to the c-axis. 

The six wires spot- welded to the sample allowed multi- 
ple redundant measurements to be performed. This im- 
proved reliability and enabled us to verify the calorime- 
try measurements using several different configurations. 
The sample resistance could be measured by a four-point 
method; a knowledge of the sample dimensions then en- 
abled the absolute resistivity to be determined to within 
10%. 

The two thermocouple junctions were formed from an 
Au/AuFe pair at either end of the sample. An alternat- 
ing resistive heating current was passed through one (to 
avoid passing the current through the sample), while the 
signal from the other was measured using a lock-in am- 
plifier. The resulting temperature oscillations serve as a 
sensitive measure of the sample heat capacity^ A sim- 
ple model of the a.c.-calorimetry system predicts the am- 
plitude and phase of the temperature oscillations (Tac) 
induced by a.c. heating: 



Tac = 



P 



K + iujC 



(3) 



where Pq is the heating power, K the thermal conduc- 
tance to the bath, C the heat capacity, and ui/2w the 
excitation frequency, assumed to be low enough that the 
thermometer can follow the temperature oscillations (the 
factor of 2ir may be assumed implicitly from this point). 
The signal therefore contains a contribution from the spe- 
cific heat and from thermal coupling to the surroundings. 

For frequencies u) 3> u> c , where lu c is the cut-off fre- 
quency K/C, the sample contribution dominates the sig- 



nal, and I Tac I can be considered to be inversely propor- 
tional to the heat capacity (which we assume to be domi- 
nated by the sample). For to <C u> c , the signal approaches 
the d.c. limit and gives a measure of the mean elevation 
of the sample temperature over that of the bath. For in- 
termediate measuring frequencies, information from the 
phase 9, can be used to extract the specific heat: 



C 



-P sm9 
u I Tac I 



(4) 



Alternatively, one can subtract a background signal taken 
at a different frequency, with 



C 



Po 



1 



1 




where ideally w% > ui c > The sample temperature 
must also be corrected for the constant d.c. component of 
the oscillatory Joule heating. This was done by repeating 
the measurement well below the cut-off frequency, also 
providing the background signal in order to estimate Cp 
using eq. JSJ. 

The cut-off frequency uj c turned out to be very tem- 
perature dependent, varying between 200Hz at 0.5K and 
over 2kHz at 1.5K, presumably due to the thermal prop- 
erties of the surrounding material. Fortunately, while 
complicating the data analysis, the reduction in lo c at 
the lowest temperatures allows the technique to be used 
down to ^100 mK. The two estimates of Cp using eqs. Q 
and JSJ are in good agreement below ~2 K. The working 
frequency L02 was generally of the order of 10 c . 

Sources of systematic error in the result might come 
from: variation of the Au Fe thermopower under pressure; 
temperature and/or frequency dependent addenda to the 
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measured specific heat due to the pressure medium, gas- 
ket and/or anvils; or any irreversibility or first-order 
character in the transitions being observed. These po- 
tential problems will be addressed in the discussion. 



III. EXPERIMENTAL RESULTS 

We present five principal results from the sample re- 
ported in this paper, and by drawing on previous work, 
we aim to place our work in a broader context. We will 
try to highlight common features found in many samples 
of CeCu2Si2, one of the defining characteristics of which 
is its variability. 

• We present the superconducting phase diagram ob- 
tained using various criteria for T c , and compare it 
to the widely quoted phase diagram determined un- 
der hydrostatic conditions by susceptibility. 

• We examine the details of the superconducting 
transition, which provides some insight into the na- 
ture of the SC state and into the sample itself. 

• We estimate the variation of the Sommerfeld coef- 
ficient 7, with pressure, and compare it to previous 
results obtained by analysis of the upper critical 
field. 

• We report the pressure dependence of the residual 
resistivity po, and exponent n, determined by a fit 
to the normal state resistivity of p = po + AT n 
(A denoting a free exponent as opposed to the 
quadratic coefficient A). A comparison of po(P) 
between different samples reveals a scaling relation 
which can be related to the theoretical enhance- 
ment of impurity scattering. 

• We explore the deviation from the scaling relation 
A oc T^ 2 , which indicates the sharp change in f- 
electron occupation number described in the intro- 
duction. The enhancement of T c and the other re- 
sults described above are shown to occur around 
the same pressure. 

Figure^shows the superconducting phase diagram de- 
termined by both resistivity and specific heat, both on 
increasing and decreasing the pressure. Two qualitatively 
different types of behavior can be seen in the same sam- 
ple, represented by the onset and completion of the re- 
sistive transition. 

If we follow the transition onset T° nset (P\ one sees 
the sharp kinks similar to those seen in Ref. [jjl (dashed 
line), along with a linear decrease of T c between 3.3 and 
4.8 GPa at a rate of 0.14 K GPa -1 . Superconductivity is 
observed however over a much smaller pressure range in 
our sample than in Ref. Il3l 

The temperature T^ =0 (P), at which the resistance 
vanishes, behaves differently from T° nset (P). It has a 
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FIG. 1: T C (P) in CeCu2Si2 determined from resistivity and 
specific heat measurements. The triangles show T c deter- 
mined from the onset of the resistive transition (T° nset ), the 
squares show its completion (T^°), and the filled circles show 
the midpoint of the specific heat jump. The numbers indicate 
the sequence of pressures. The dotted line shows T c deter- 
mined by susceptibility in a different sample, also in a helium 
pressure mediums- 



narrower peak with a maximum at slightly higher pres- 
sure. T c fl=0 agrees closely however with the transition 
seen in the specific heat (see below). When a magnetic 
field was applied, T c fl=0 and the specific heat anomaly 
shifted in agreement. 

The large resistive transition widths found in 
CeCu2Si2 at high pressure are often blamed on a lack 
of hydrostaticity due to the pressure medium. As helium 
was used in this case, we can rule out pressure inhomo- 
geneities and concentrate on the sample itself. Further 
information about the SC state comes from the effect of 
measurement current on the transition width. For ex- 
ample, at 1.78 GPa high current led to the upper part 
of the transition disappearing, and a resistive transition 
can even be recovered with a narrow width comparable to 
that close to ambient pressure. This is presumably due 
to the presence of filamentary superconductivity, with 
a higher T c , whose critical current density is exceeded. 
These broad resistive transitions appear to be a universal 
feature of CeCu2Si2 at high pressure. Let us recall that 
even for the highest T° nset measured in a single crystal, at 
2.4K, a tail of 1% of the normal state resistivity remained 
well below 2 K, vanishing only at 1.5 K.H The status of 
the superconductivity of CeCu2Si2 between T° nsct and 
T^- =0 remains mysterious. 

Figure [3 compares the superconducting transition in 
resistivity and specific heat at three different pressures. 
At 2.38 GPa the resistive transition is broad and the 
sharp specific heat jump at 0.73 K begins at the point 
where the resistance falls to zero. At 3.67GPa the specific 
heat jump, at 1.35 K, is much larger, and remains sharp 
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FIG. 2: Superconducting transition at three pressures in (a) 
resistivity and (b) specific heat. Note the width of the re- 
sistive transitions, and the fact that the start of the jump 
in specific heat coincides with the completion of the resistive 
transition. 



(and did so at intervening pressures), while the corre- 
sponding resistive transition has narrowed considerably. 
At 4.07 GPa (not shown) where T*=° has a maximum 
around 1.6 K in both p and Cp, the specific heat peak 
has already started to broaden and collapse in amplitude, 
while at the same pressure the resistive transition is at its 
narrowest since ambient pressure. As T c is driven to zero 
at high pressure, the superconducting Cp jump becomes 
smaller and broader (as shown at 4.67 GPa) until it is no 
longer visible. When the pressure was reduced, the Cp 
peak recovered its shape, indicating the reversibility of 
the bulk pressure-induced behavior. 

The dramatic increase in the apparent size of the 
superconducting jump is intriguing, and might suggest 
the presence of strong coupling, 35 or other qualitative 
change in the SC state. Although the apparent value of 
(ACp/7T) Tc is clearly less than the BCS ratio of 1.43, 
similar a.c. measurements on CeCoIns in an argon pres- 
sure medium indicate that there is a substantial contri- 
bution to the measured heat capacity from addenda^ In 
helium we would expect this to be even more significant 

The increase in the Cp jump size might itself be 
an artefact of the uncalibratcd a.c.-calorimetry method; 
nevertheless (ACp /jT) T does appear to show a max- 
imum at a pressure coinciding with the increase in T c . 
Furthermore, the assumption of strong coupling provided 
the best fit to H c % for measurements of the upper critical 
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FIG. 3: Estimate 'y(P) of the Sommerfeld coefficient from a.c.- 
calorimetry signal at 2 K (triangles), scaled for comparison 
with that deduced from H C 2 measurements (squares)^- The 
lines are guides for the eye. 



field in another sampled 

The electronic specific heat coefficient 7, and hence 
the effective mass m*/m, can be estimated by following 
the calorimetric signal C/T, at a fixed temperature and 
measurement frequency above the superconducting tran- 
sition, though this includes constant or slowly varying 
addenda from the helium, diamonds etc. Figure y] shows 
the estimate j(P), along with the value deduced from 
measurements of the upper critical field in Ref. ^3 A 
single constant scale factor has been introduced, show- 
ing that the two curves can be superimposed. There is a 
clear anomaly in 7 at 4 GPa (just below the pressure cor- 
responding to T™ ax ), superimposed on a constant reduc- 
tion with pressure. The effective mass is also reflected in 
the initial slope of the upper critical field H' c2 (T c ), which 
in our sample also had a maximum at the same pressure 
as the peak in 7. 

The residual resistivity po, has a huge peak at a pres- 
sure slightly higher than the maximum in T c . The mag- 
nitude of this peak varies by a factor of more than ten 
between samplesi^Si However, it is possible to scale the 
residual resistivities from different samples onto the same 
lorentzian curve (see Fig-HJ. A constant value p$, differ- 
ent for each sample, is subtracted from po at each pres- 
sure, and the result is multiplied by a scale factor, a (i.e. 
Po = a(po — po)), so that all lie on the curve defined by 
sample SI, which has the highest residual resistivity (i.e. 
asi = I). 

According to the theoretical prediction, the residual 
resistivity po is given as2i 



po = Bn imp \u(0)\ 2 In 



drii 
def 



unit 

Po ) 



(6) 



where the coefficient B depends on the band structure of 
host metals, n.i mp is the concentration of impurities with 
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FIG. 4: Enhancement of residual resistivity in several differ- 
ent CeCii2Si2 samples, scaled to a universal pressure depen- 
dence, with a and pg being normalizing factors. The maxi- 
mum in po is at a pressure slightly higher than that corre- 
sponding to the maximum in T c . The inset shows n for a fit 
of p — po + AT n . Filled squares indicate results from this 
work, the rest are reported in Ref. [23 (S and C refer to the 
original labels and a retained for continuity). 



moderate scattering potential u(q) coming from disorder 
other than Ce ions, iVp is the density of states of quasi- 
particles around the Fermi level, and the last term rep- 
resents the residual resistivity due to unitary scattering 
mainly arising from any deficit or defect of the Ce ions. 
The scaling behavior of po shown in Fig. 0] would be possi- 
ble if the universal form is given by In \ (—drif / dei) ^/N-f]. 
It is an open question whether the observed lorentzian 
form is indeed reproduced by the theory of Ref. |U 

There is a striking correlation between the scaling fac- 
tor a and the behavior of T c . The sample measured in 
helium reported in this paper and sample CI, pressur- 
ized in steatite, both have similar values of a, and both 
yonset an( j jiR=o a g ree over almost the entire pressure 

range. Sample SI, with the highest po at P v , has a lower 
ymax ^ -yfi K), and the superconductivity disappears 
at a lower pressure. Samples C2 and S2 have scaling 
factors a around 14, and show a higher maximum T c , 
with superconductivity extended over a greater pressure 
range than in the samples with larger residual resistivi- 
ties. These differences between samples, both in po and 
T c , are vastly amplified from their appearance at ambient 
pressure. 

The inset in Fig. 0] shows the result of a fit to p = 
Po + AT n above T c . There are two important points to 
note here. Firstly, at the pressure slightly higher than the 
maximum T c , p(T) is linear in T up to about 25 K. Sec- 
ondly, the exponent appears surprisingly large (n ~ 2.7) 
at the slightly higher pressure corresponding to the max- 
imum po- This is difficult to understand without taking 
into account the resistivity due to impurity scattering. 
In sample SI, reported in Ref. l26t the residual resistivity 
reaches ~ 160 pQcm at P v , compared to a maximum 



of 35 pflcm for the sample reported here. p(T) then 
showed a falloff with temperature very similar to that of a 
Kondo impurity system. In other samples, this behavior 
is hidden by the usual positive temperature dependence 
of the resistivity. Contrary to the usual situation, where 
the lowest po possible is sought, this example shows how 
samples whose residual resistivities are large at ambient 
pressure can reveal interesting physics at high pressure. 
Even if a negative temperature dependence is not seen, 
the power-law fit to the resistivity is affected, deviat- 
ing from the linear relationship predicted in section IIVI 
and leading to anomalous values of n. At lower pres- 
sure, the A coefficient is an order of magnitude larger, so 
(for example) almost linear resistivity is observed at the 
pressure corresponding to T c max . Note that a quadratic 
temperature dependence of p was recovered at the lowest 
temperatures when superconductivity was suppressed by 
a magnetic field greater that H C 2- 

The normal state resistivity of heavy fermions can usu- 
ally be understood in terms of the Kondo lattice models 
At high temperature the f-electron moments are localized 
and disordered, the resistivity is large and dominated by 
the scattering from spin disorder, with a characteristic 
— InT slope. As the temperature is reduced, Kondo sin- 
glets form below a characteristic temperature Tk, and 
coherence effects in a periodic lattice cause the resistivity 
to drop below a maximum, at T™^, which can be con- 
sidered as proportional to Tk- For T <C Tk away from 
a critical point, Fermi-liquid-like behavior is recovered, 
with p ~ AT 2 , where A oc T^ 2 and reflects the hugely en- 
hanced effective mass caused by interactions between the 
f-electrons. In a real system where Tk is not too large, 
a second peak in the resistivity occurs at T™ ax > T™, 
due to the crystalline electric field (CEF) effect 38 i 39 (see 
inset of Fig. The low temperature behavior then re- 
flects the characteristics of the lowest CEF-split f-level. 
When pressure is applied, T™ ax remains fairly constant, 
while Tk rapidly increases, seen via the rise in 7 1 1 max . 
When T K > A CE f (Aqef is the CEF splitting between 
the ground and excited states) the full 6-fold degener- 
acy of the J=5/2 4f 1 multiplet is recovered, even at the 
lowest temperatures. As a result the resistivity maxima 
at T™ ax and T™ ax merge into a single peakJ-S Similar 
behavior in the magnetic component of the resistivity is 
found in all Ce compounds studied (such as CeCusAurfl 
CePd 2 Si 2! 22 CePdaGe^). 

In Fig. [5] the A vs T™ ax scaling is explored in both 
CeCu2Si2 and CeCu2Ge2- The value of A was deter- 
mined from the slope of the normal state resistivity ver- 
sus T 2 , despite the non- Fermi liquid behavior shown in 
the inset of Fig.^J However, if one allows the exponent n 
to vary between 1 and 2, the resulting coefficient will not 
vary more than a factor of two, which is within the scatter 
of the data. There are two regions where the predicted 
A oc (7« lax )- 2 relationship is followed, separated by an 
abrupt drop in A of over an order of magnitude. The 
collapse of A seems closely connected with the enhance- 
ment of superconductivity, it is at the start of this drop 
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FIG. 5: Plotted against T™ ax (defined in inset), a measure of 
the characteristic energy scale of the system, are (a) the bulk 
superconducting transition temperature, (b) the residual re- 
sistivity and estimate 7 of the Sommerfeld coefficient, and (c) 
the coefficient A of the p ~ AT 2 law of resistivity, including 
data from CeCu2Ge2. Note the straight lines where the ex- 
pected A oc (7\ max )- 2 scaling is followed. The maximum of 
T c coincides with the start of the region where the scaling re- 
lation is broken, while the maximum in residual resistivity is 
situated in the middle of the collapse in A. Pressure increases 
towards the right-hand side of the scale (high T 1 ™ 5 " 1 ). 



that T c has maximum, and the superconductivity has dis- 
appeared by the point where the A oc (y™ ax ) -2 scaling 
is recovered. The residual resistivity however, peaks at 
around the midpoint of the drop in A, and this is the 
point where P v is defined. 



IV. THEORY OF T-LINEAR RESISTIVITY AND 
ENHANCED SOMMERFELD COEFFICIENT 

Various unconventional properties observed around 
P ~ P v have been explained, at least qualitatively, by 
a series of theoretical investigations on the basis of an 
extended Anderson lattice model However, the T- 
linear temperature dependence of the resistivity observed 



in a narrow region around P ~ P v remains as yet unex- 
plained. In Ref. l22t microscopic calculations showed that 
the static limit of the effective interaction r^°'(q) between 
quasiparticles is enhanced greatly around P ~ P v , and is 
almost independent of q, the momentum transfer, up to 
~3/2 of pf, reflecting the local nature of critical valence 
fluctuations. This implies that the valence fluctuation re- 
sponse function Xv(q,w), is also almost ^-independent in 
the low frequency region. Based on this observation, we 
present here a phcnomcnological theory explaining the T- 
linear resistivity and the enhancement of the Sommerfeld 
coefficient 7 around P ~ P v . 

We adopt an exponentially decaying phcnomcnological 
form for the valence-fluctuation propagator (dynamical 
valence susceptibility) Xv^ 



Xv(g,w) = i / dte ,ut <[n f (( Z ,t),Tif(-g ) 0)]) (7) 
Jo 



K 



— , for q < q c ~ pf 



(8) 



where rif(g) is the Fourier component of the number 
of f-electron per Ce site, if is a constant of 0(1), 
and tu v parameterizes the closeness to criticality. lo v 
is inversely proportional to the valence susceptibility 
Xv (0,0) = -(dm/defy. 

The real and imaginary parts of the retarded self- 
energy Syf(p, e + ifl) respectively give a measure of the 
quasiparticle effective mass and lifetime. They can be 
calculated using a simple one-fluctuation mode exchange 
process (see Fig. [SJ and given as follows: 



Re££(p,e) 



•— y iai 2 



dx- 



coth 



2T 



tanh^a 



, (9) 



p-q 



ImE*(p,e) 



K 



p-q 



^ + (e - e 



p-q 



x coth • 



p-q 



2T 



tanh^za 
2T 



■ (10) 



where A is the coupling between quasiparticles and the 
valence fluctuation modes, and £ p is the dispersion of the 
quasiparticle. For simplicity, A is assumed to be constant 
without wavenumber or frequency dependence. 

In typical limiting cases, I jlOjl can be straightforwardly 
calculated in the approximation £p-q — — vq cos 9, where 
is the angle between p and q, v is the quasiparticle 
velocity, and p is assumed to be on the Fermi surface, 
i.e., p = pf- 



T = 0, e ± 0: 

ImS^ f (p F ,e) 



\M 2 Kq 2 c 
32n 2 v 



In 1 + ^ 



(11) 
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p-q 



FIG. 6: Feynman diagram for the self-energy given by eqs. @ 
and HUH . The solid line represents the Green function of 
the quasiparticles, the wavy line the propagator of the va- 
lence fluctuations, and the filled circle the coupling between 
valence-fluctuation modes and the quasiparticles. e n and ui m 
are the Matsubara frequency of the quasiparticle and fluctu- 
ation propagators, respectively. 



where q c is the cutoff wavenumber of the order of fc F . 



e = 0, < T< e F : 
Im£* f (p F ,0) - 



\M 2 k 

8tt 2 v 

y 



dqq 



vq/2T 



vq/2T 



dy x 



-(cothf -tanh|)(12) 
(uj v /Ty + y 2 2 2 

where y = vqcos9/2T. Since vq 3> T holds in the domi- 
nant region of q-space, the integration with respect to y 
can be performed, to a good accuracy, leading to 



Im£? f (p F) 0) 



' A ' 2 ^ T tan^, (13) 



4ir 2 v u> v 



where we have made approximation that the range of 
integration is restricted as — 1 < y < 1 in which the last 
factor in 1|12|) is approximated as 2/y. Then, 



Im£S(pF,0) 



|A| 2 ifg c 2 

Al7 2 V 




T< LU V 
T>CJ V 



(14) 



The latter result, Im£ v f (p F , e = 0) oc T/lo v for T ^ u> v , 
implies that almost all the critical valence-fluctuation 
modes can be regarded as classical at T > w v , and T- 
linear dependence stems from the asymptotic form of 
coth(a;/2T) ~ 2T/x, essentially the classical approxima- 
tion of the Bose distribution function. 

The real part of the self-energy, JSJl, can be calculated 
easily at T — and e ~ 0, leading to 

ReE* f (p F ,e) - S^(p F ,0) 



W 2 Ke 
4tt 2 



(vqt)'- 



dqq 2 



■In 



dt 



vqt 



2uol 



(vqty 



■In 



TTU) v vq\t\ 

R 



vqt 



(vqt) 2 } 2 
, (15) 



where t — cos 9. In the limit lu v <C vpp, integration with 
respect to t in (|15fl leads to 



Re£*(p F , £ ) 



\X\ 2 Kq 2 e f 1 1-u 2 

' du— In 



2n 2 v lo- 



v Jo 



0J V 



-(16) 
u 

(17) 



where u = vqt/u> v . 

The T-linear dependence of Im£^ f (p, 0), for T > ui v , 
(|14|l . implies T-linear resistivity, as the quasiparticles 
are subject to the large angle scattering by the criti- 
cal valence-fluctuation modes. These are effective in a 
wide region in the Brillouin zone due to their local nature 
and easily couple to the Umklapp process of quasiparticle 
scattering. This result is consistent with the experimen- 
tal fact that T-linear resistivity is observed in a narrow 
pressure region around P v , which is considered to cor- 
respond to a nearly critical valence transition of the Ce 
ion. 

Such a T-linear dependence has been discussed in the 
context of high-T c cuprates with a marginal Fermi liq- 
uid (MFL) assumption^ and charge transfer fluctua- 
tions were once considered as an origin for MFL 
while further theoretical models have been put forth up 
to nowm^ Excepting the T-linear resistivity, the present 
result is different from MFL behavior. The self-energy 
exhibits different energy dependence, while the idea for 
the origin of our singular behavior shares aspects simi- 
lar to the first idea of a charge transfer mechanism for 
high-T c cupratesi^ii £(e) in the MFL model is given 
as E(e) oc (elne — i|e|)^ which is indeed different from 
the present case [Eas. <|ll|l and (|17fl ]. In any case, it is 
to be noted that T-linear resistivity is accompanied by 
the peak of T c in both systems, high-T c cuprates and 
CeCu 2 Si 2 . 

The result (|17fl implies that the mass enhancement ( 1 — 
9ReS^ f (e) / de) is expected around P ~ P v . Namely, the 
effective mass is given by 



m oc m- 



(18) 



where fh is the effective mass renormalized by the conven- 
tional correlation effect, leading to heavy electrons, i.e. 
not including the effect of critical valence fluctuations. 
This latter effective mass m exhibits a drastic decrease 
around P ~ P v , while the second factor in i|18fl is en- 
hanced. Both effects should be reflected in the Sommer- 
feld coefficient 7, so that the peak of 7 oc m* is shifted to 
the lower pressure (larger fh) side, and the anomaly of 7 
due to the valence fluctuations may be smeared to some 
extent. Nevertheless, some trace should be observed. 
(The shift of peak of 7 can be understood as the super- 
position of the two trends using a model P-dependcncc 
of fh and lo v .) Indeed, the present experimental result 
presented in Figs. and |S] may be explained by this ef- 
fect. 
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V. DISCUSSION 

Our calorimetric results in such extreme conditions de- 
serve some discussion, in particular the considerable ap- 
parent increase in the specific heat jump at the supercon- 
ducting transition when P v is approached. A very large 
specific heat jump at T c would be strongly reminiscent 
of the huge value found in CeCoIxis. 46 It is therefore a 
legitimate question to ask how much the results of the 
uncalibrated AC calorimetry technique under pressure 
can be relied on to give an accurate measurement of the 
specific heat. 

The model used to extract the specific heat from the 
amplitude and phase of the temperature oscillations takes 
no account of the heat capacity of the solid helium, di- 
amonds, or surrounding pressure apparatus, or the es- 
sentially three-dimensional nature of the situation. Sec- 
ondly, the thermopower of the Au Fe thermocouple has 
been assumed not to vary with pressure (Ref. 40 indi- 
cates that it varies by no more than 20% up to 12 GPa). 

Nevertheless, the superconducting transition observed 
corresponds to ~ 100% of the signal amplitude, indicat- 
ing that the addenda are a minority contribution to the 
total signal. Runs at several different frequencies agree to 
within 10-20% after the amplitude and phase are com- 
bined, with the discrepancy possibly due to frequency- 
dependent addenda. Kapitza resistance between the 
sample and helium is likely to better decouple the sam- 
ple from its surroundings at very low temperature. If 
the specific heat is calculated using the two-frequency 
method [eq. J5J], the result agrees (< 5%) with that cal- 
culated using the amplitude and phase up to at least 2T C . 
Given these observations, it seems reasonable to accept 
our results as a good first approximation to Cp, to within 
a constant scaling factor, and with an unknown but rel- 
atively small component due to addenda. 

Furthermore, the apparent anomaly in the normal 
state specific heat shown in Fig.|21was measured at a fixed 
temperature and frequency above the superconducting 
transition, with pressure the only independent variable. 
The small peak in 7 is consistent with the maximum in 
the initial slope of the upper critical field observed at the 
same pressure, though the interpretation of the latter de- 
pends on whether the sample can be considered to be in 
the clean or dirty limit, or somewhere in between. 

Having addressed the experimental questions, let us 
discuss some other remaining points. The merging of 
rpmax anc j ymax seems ^ De De a general feature at P v in 

compounds where a critical valence transition is thought 
to exist. It can be understood as follows: 

The so-called Kondo temperature Tk, related to Xjj max 
(i = 1,2), depends crucially on the degeneracy (2£ + 
1) of the local f-state: T K ~ Dexp[-l/(2l + 1)pf\J\], 
where D is the bandwidth of conduction electrons, /?f 
the density of states of conduction electrons at the Fermi 
level, and J the c-f exchange coupling constant^ Even 
though the 6-fold degeneracy of the 4f-state is lifted by 
the CEF effect, leaving the Kramers doublet ground state 



and excited CEF levels with excitation energy AceFj the 
Kondo temperature Tk is still enhanced considerably by 
the effect of the excited CEF levels^ 

The technical degeneracy relevant to the Kondo effect 
is affected by the broadening AE of the lowest CEF level. 
If AE <C Acef, the degeneracy relevant to Tk is 2-fold. 
On the other hand, if AE > Acef, it increases to 4- or 
6-fold. The level broadening is given by AE ~ zitpy\V\ 2 
where \ V\ is the strength of c-f hybridization, and z is the 
renormalization factor which gives the inverse of mass 
enhancement in the case of a lattice system. It is crucial 
that AE is very sensitive to the valence of Ce ion because 
z is essentially given by q [eq.Q]. In particular, the 
factor z increases from a tiny value in the Kondo regime, 
2 ~ (1 — Bf) < 1, and approaches unity in the so-called 
valence fluctuation regime. 

Since the factor | ^| 2 3> Acef in general for Ce- 
based heavy electron systems, the ratio AE'/Acef, which 
is much smaller than 1 in the Kondo regime, greatly ex- 
ceeds 1 across the valence transformation around P ~ P v , 
leading to the increase of the technical degeneracy of f- 
state, irrespective of the sharpness of the valence transfor- 
mation. Therefore, T" lax should merge with T™^, which 
corresponds to 4- or 6-fold degeneracy of 4f-state due to 
the effect of finite temperature, i.e., T ~ Acef- This 
may be the reason why T™ ax increases and approaches 
ymax p ressure where T c exhibits the maximum, and 
the KW ratio changes between strongly and weakly cor- 
related classes. 

While the experimental picture of CeCu2Si2 presented 
in this paper is more complete than the theoretical, a 
large number of the features found around P v follow di- 
rectly from the valence fluctuation approach and the ad- 
dition of a U c f term to the hamiltonian. The linear resis- 
tivity is explained in section llVl as is the local maximum 
in the electronic specific heat, possibly due to the renor- 
malisation of the effective mass due to valence fluctua- 
tions, superimposed on an overall decrease with pressure. 
The enhancement of the residual resistivity at low tem- 
perature follows from the renormalisation of impurity po- 
tentials by valence fluctuations. The relative positions of 
the peaks in T c , 7, and po are consistent with the valence 
fluctuation scenario, but for a more precise comparison 
more detailed calculation would be needed. 

Other features yet to be fully addressed with the cur- 
rent model are observed to occur in the valence fluctua- 
tion region. They are the apparent increase in the specific 
heat jump at T c , the temperature dependence of the im- 
purity contribution to the resistivity, and the nature of 
the superconducting state between the onset and com- 
pletion of the superconducting transition. 

The presence, and indeed enhancement, of supercon- 
ductivity so far from the disappearance of magnetic order 
calls into question whether magnetic mediation is really 
the sole mechanism of superconductivity in CeCu2Si2- 
The evidence presented here, along with other anoma- 
lous behavior seen at a pressure well separated from the 
disappearance of magnetism, strongly suggests the pres- 
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FIG. 7: Schematic P-T phase diagram for CeCu2(Si/Ge)2 
showing the two critical pressures P c and P v . At P c , where 
the antiferromagnetic ordering temperature Tjv — > 0, super- 
conductivity in region SC I is mediated by antiferromagnetic 
spin fluctuations; around P v , in the region SC II, valence 
fluctuations provide the pairing mechanism and the resistiv- 
ity is linear in temperature. The Temperatures y 1 max i an d 
T™ ax , merge at a pressure coinciding with P v . 



ence of a second quantum critical point in CeCu2Si2, this 
time related to quantum fluctuations between electronic 
configurations rather than to collective spin instabilities. 
While magnetic pairing may be responsible for supercon- 
ductivity at the magnetic QCP, critical valence fluctua- 
tions are responsible for pairing at P v . The recent result 
in CeCu2(Sio.gGeo.i)2r^ where two separate peaks of T c 
are observed, suggests the validity of the present point of 
view. 

Figure [7| shows a schematic phase diagram for the 
CeCu2(Si/Ge)2 system. The two critical pressures, P c 
and P v , are respectively defined by the disappearance 
of magnetic order as T/v tends to zero, and by the re- 
gion of linear resistivity where po has a maximum and 
jimax ymax^ accom p an i ec [ by a maximum in T c . In 

CeCu2Si2 and CeCu2Ge2 the two critical pressures are 
widely separated. In compounds such as CePd2Si2 on 
the other hand, superconductivity is found in a narrow 
pocket, seemingly directly connected to the disappear- 
ance of magnetism as TV — > 0. However, many of the 
other anomalies listed in tableware still observed in this 
system, and are difficult to explain within a purely spin 
fluctuation picture. If a valence instability is present in 
CePd2Si2, P v is superimposed on P c , as identified by the 
pressure at which T™ ax ~ T 2 max i2S The physics associated 
with valence change in CeCu2Si2 may thus also play an 
important role in other heavy fermion superconductors. 



Linear resistivity and an enhancement of po have also 
been seen in the CeTIns compounds^^ where T is Co, 
Rh or Ir. For this family, superconductivity extends over 
a relatively broad pressure range, and it may be that va- 
lence fluctuations also play a role with a critical valence 
pressure separate from any magnetic instability. 

Finally, it is worth addressing the physical interpre- 
tation of the valence fluctuation mediated pairing inter- 
action. We emphasize that this intuitive explanation is 
rather speculative, but we think that it is sufficiently use- 
ful to merit inclusion. 

A clue comes from the likely nearest neighbor pairing, 
implied by the largely local nature of the interaction, 
and the prediction of d-wave pairing symmetry. One can 
imagine an almost filled f-band, with each occupied f 1 site 
experiencing a Coulomb repulsion U c f from the respective 
conduction electrons. As the pressure is increased and ef 
moves closer to the Fermi level eF , there will come a point 
where ef + U c f — e-p and the f-band will start to empty. 
On an individual 4f° 'hole' site, the U c i interaction will 
be absent, thus an increased density of conduction elec- 
trons would be energetically favorable at this position. 
If this extra 'screening' conduction electron density is 
not strictly localized onto the atom itself, but spills onto 
neighboring sites, the f-electrons on Ce atoms around the 
original 'hole' site will feel an increased repulsion. The 
tendency to transfer electrons from the f to conduction 
bands will be locally reinforced, explaining intuitively the 
increasingly catastrophic drop in nf for larger U c {, pre- 
dicted in Ref.|23 For large enough U c f, phase separation 
would be expected to occur for some values of ef. 

The attractive pairing interaction can be understood as 
follows: Consider an isolated pair of 4f° 'holes', accom- 
panied by their cloud of conduction electrons. If these 
are separated by two lattice positions, with an interven- 
ing filled 4f x site, the two clouds of conduction electrons 
will overlap at the intermediate site, further increasing 
the Coulomb energy at that point. It would therefore 
be energetically favorable for the two 'holes' to be on 
neighboring atoms, thus the attractive interaction. The 
attractive interaction between 'holes' is equivalent to that 
between 'electrons', so that this argument would give 
an intuitive understanding of the origin of the valence- 
fluctuation mechanism of superconductivity. 



VI. CONCLUSIONS 

The enhancement of superconductivity in 
CeCu2Si2 under pressure is found to coincide with 
a number of anomalies in the superconducting and 
normal state properties that are hard to explain in 
a purely spin fluctuation scenario. Many of these 
anomalies are directly related to an abrupt change in 
valence of the Ce ion, while others can be indirectly 
connected to such a transition. We propose a second 
critical pressure P v at around 4.5 GPa where critical 
valence fluctuations provide the superconducting pairing 
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mechanism. An extended Anderson lattice model 
with Coulomb repulsion between the conduction and 
f-electrons predicts an abrupt change in Ce f-level 
occupation. The associated fluctuations are sufficient to 
explain the observed enhancement of T c , the T-lincar 
normal state resistivity, the enhancement of the residual 
resistivity, and the peak in the electronic specific heat 
coefficient 7. 
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